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I discuss the applicability of variational methods to the study of non-pcrturbativc 
aspects of QCD. An illustration of the capabilities of the method pioneered by Ko- 
gan and Kovner is given through the analysis of the deconfinement phase transition 
in gluodynamics in 3+1 dimensions. 



1. Introduction 

There is very little doubt that Quantum Chromodynamics (QCD) is the 
correct theory of the strong interactions. Understanding the properties of 
the vacuum sector of QCD - phenomena such as confinement and chiral 
symmetry breaking - remains, in spite of years of attempts, one of the 
main problems in modern quantum field theory. 

Amongst the variety of approaches used to tackle these problems, there 
have been several attempts to apply versions of the variational Rayleigh- 
Ritz method [1]. The crux of my argument is that such methods, in partic- 
ular Kogan and Kovner's version [1,2], offer a possibility to analytically 
obtain non-perturbative dynamical information directly from QCD. Al- 
though this information is partial and incomplete, the variational approach 
provides significant qualitative understanding, and, to a limited extend, 
quantitative estimates. 

When attempting to apply the variational principle to quantum field 
theory (QFT) one is faced with the discussed by Feynman [3]. In what 
follows I will describe the construction of a variational trial state which, 
to a certain extent, overcomes these difficulties, and then illustrate the 
variational approach's capabilities by discussing its application to the study 
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of the dcconfining phase transition in gluodynamics in 3+1 dimensions [4,5]. 
2. The variational Ansatz 

The first criterion to be considered when constructing a trial variational 
state is that the trial state ought to be general enough to allow, through 
variation of its parameters, the relevant physics to be spanned. 

Then there is the problem of calculability. In QFT, our ability to eval- 
uate path integrals is, to say the least, limited. The requirement of calcu- 
lability almost unavoidably restricts the possible form of the trial state to 
a Gaussian. 

Another serious problem is that of "ultraviolet modes" . The main mo- 
tivation of a variational calculation in a strongly interacting theory is to 
learn about the distribution of the low momentum modes of the field in the 
vacuum wave functional. However, the VEV of the energy (and all other 
intensive quantities) is dominated by contributions of high momentum fluc- 
tuations, for the simple reason that there are infinitely more ultraviolet 
modes than modes with low momentum. Therefore, even if one has a clear 
idea of how the WF at low momenta should look, if the ultraviolet part 
of the trial state is even slightly incorrect the minimization of energy may 
lead to absurd results. Due to the interaction between the high and low 
momentum modes, there is a good chance that the infrared (IR) variational 
parameters will be driven to values which minimize the interaction energy, 
and have nothing to do with the dynamics of the low momentum modes 
themselves. 

Finally, in gauge theories there is the additional complication of gauge 
invariance. Allowable wave functional must be invariant under the time 
independent gauge transformations. If one does not impose the Gauss' 
law on the states exactly, one is not solving the right problem. The QCD 
Hamiltonian is only defined on the gauge invariant states, and its action on 
non gauge invariant states can be modified at will. Thus, by minimizing 
a particularly chosen Hamiltonian without properly restricting the set of 
allowed states, one is taking the risk of finding a "vacuum" which has 
nothing to do with the physical one, but is only picked due to a specific 
form of the action of the Hamiltonian outside the physical subspace. 

The trial variational state proposed by Kogan and Kovner [2] explicitly 
satisfies Gauss' law, has the correct UV behaviour built-in, and allows for 
analytical calculations - yielding non-perturbative results - to be carried 
out. 
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In order to illustrate the capabilities of this variational method, I now 
turn to the analysis of the deconfinent transition in gluodynamics. 

3. Deconfinement phase transition 

In the almost 25 years since the pioneering work of Polyakov [6] and 
Susskind [7] much effort has been devoted to attempts to understand both 
the basic physics and quantitave features of the deconfining phase transition 
of QCD. 

The high temperature phase is becoming well understood, and is widely 
believed to resemble a plasma of almost free quarks and gluons. However, 
the transition region, T c < T < 2T C , is very poorly understood. This 
region is the most interesting since it is there that the transition between 
'hadronic' and 'partonic' degrees of freedom occurs. 

The study of the transition region is a complicated and inherently non- 
perturbative problem which has mostly evaded treatment by analytical 
methods. Recently, the method introduced by Kogan and Kovner [2] was 
applied to the study of the phase transition in SU(N) gluodynamics [4,5]. 
Following [4,5], we minimize the relevant thermodynamic potential at finite 
temperature, i.e. the Helmholtz free energy, on a set of suitably chosen trial 
density matrix functionals. We consider density matrices which, in the field 
basis, have Gaussian matrix elements and where gauge invariance is explic- 
itly imposed by projection onto the gauge-invariant sector of the Hilbert 
space 

p[A,A'\ = j DU expj - ± jf yi{x)Gr^\x,y)A){y) 

+ A' i Ua (x)GT j lab (x,y)A' J ub (y) 2A?(x)H? b (x,y)A' 3 ub ( y )] J , (1) 

where f xy = Jd 3 xd 3 y 7 DU is the SU(N) group-invariant measure, and 
under an SU (N) gauge transformation U 

A${x) - AY a {x) = S ab (x)A b (x) + \«{x) , (2) 

with 

S ab = itr(r a C/W) , A" = -tT{ T a U^diU) , (3) 
2 9 

and ^ f° rm an iV x iV Hermitian representation of SU(N): [-rr,-y] = 
if^cil with normalization tr(r a r b ) = 2S ab . 
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We take the variational functions diagonal in both colour and Lorentz 
indices, and translationally invariant 

G-^ix, y) = S ab S i jG~ 1 (x - y) , H$(x, y) = 5 ab S l3 H(x - y) . (4) 

Further, we split the momenta into high and low modes with k ^ M and 
restrict the kernels G _1 and H to the one parameter momentum space 
forms 



H, k < M 

0, k > M ■ W 



The logic behind this choice of ansatz is the following. At finite temper- 
ature we expect H(k) to be roughly proportional to the Bolzmann factor 
cxp{ — E(fc)/3}. In our ansatz, the role of one particle energy is played by the 
variational function G~ l {k) and its form is motivated by the propagator 
of a massive scalar field, i.e. (fc 2 + M 2 ) 1 / 2 . We will be interested only in 
temperatures close to the phase transition, and those we anticipate to be 
small, T c < M. For those temperatures one particle modes with momenta 
k > M are not populated, and we thus put H{k) = 0. For k < M the 
Bolzmann factor is non- vanishing, but small. Further, it depends only very 
weakly on the value of the momentum. With the above restrictions on the 
kernels, only two variational parameters, M and H, remain. 

Importantly, the density matrix functional in eq. (1) describes, for 
H = 0, a pure state p = ^[A] >< *&[A]\ where \I>L4] are Gaussian wave 
functionals. For H ^ 0, eq. (1) describes a mixed state with \H\ propor- 
tional to the entropy of the trial density matrix. 

The expectation value of a gauge invariant operator in the variational 
state eq. (1) is then given by 

{0) A ,u - Z^TripO) 

= Z- 1 J VUVAO(A,A')- 

• exp | - X - \aG~ x A + A^G- 1 A lU - 2AHA' U ] J 

where Z is the normalization of the trial density matrix p, i.e. 

Z = Tip = J VUVA exp j - ^ AG~ 1 A + A U G~ 1 A U - 2AHA U 

To evaluate the above expressions we first perform, for fixed U(x), the 
Gaussian integration over the vector potential A. For Z we get, in leading 



, (6) 

A'=A 



(7) 
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order in H , 

Z = |wcxp|-iA(^ + |(5 + ^))A+^Gtr(5 + 5 T )}. (8) 

We now integrate out the high momentum, k 2 > M 2 , modes of U per- 
turbatively to one-loop order. This effects a renormalization group trans- 
formation on the low modes, replacing the bare coupling g 2 by the running 
coupling g 2 (M) [2,9]. To one- loop accuracy, the coupling g 2 (M) runs iden- 
tically to the Yang-Mills coupling [9] . 

The normalization Z can be then interpreted as the generating func- 
tional 

Z = Trp = J VUe- s ^ (9) 

for an effective non-linear cr-model for the low momentum modes (k 2 < M 2 ) 
defined by the action 

{U ] dU - dU^U){dUU^ - UdU^) 

- -L-HMhrUhrU , (10) 

where U independent pieces have been dropped. 

The matrix U plays the same role as Polyakov's loop P at finite tem- 
perature — the functional integration over U projects out the physical sub- 
space of the large Hilbert space on which the Hamiltonian of gluodynamics 
is defined. 

This cr-model has a phase transition [2] at the critical point (for Aqcd = 
150Mev, N = 3 and with the one-loop Yang-Mills (3 function) 

M c = A QCD eTT = 8.86A QC D = 1.33Gev. (11) 

For M < M c , the cr-model is in a disordered, SU (N) L !g>SU (N) R symmetric, 
phase with massive excitations and where (U) = 0. Since U is the Polyakov 
loop, this corresponds to a confined state. When M > M c , the cr-model is 
in an ordered, SU(N)v symmetric, phase with massless Goldstone bosons 
for which (U) ^ 0, corresponding to a deconfined state. With this analysis 
we have established a correspondence between the cr-model phase transition 
and the deconfincment transition in SU (N) gluodynamics. 

In fact, this correspondence can be argued in rather general terms. Let 
us ask ourselves what would happen if we did not restrict H to be small, 
and more generally did not restrict the functional forms of G(k) and H(k) 
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in our variational ansatz. We could still carry on our calculation for a while. 
Namely we would be able to integrate over the vector potentials in all aver- 
ages, and would reduce the calculation to a consideration of some non-linear 
cr-model of the J7-ficld. This cr-model quite generally will have a symmetry 
breaking phase transition as the variational functions G(k) and H(k) are 
varied. Since at this transition the Polyakov loop U changes its behavior, 
the disordered phase of the cr-model corresponds to the confining phase of 
the Yang-Mills theory, while the ordered phase of the cr-model represents 
the deconfincd phase. Thus, in order to study deconfinement in the SU(N) 
Yang-Mills theory, we should analyze the physics of each cr-model phase 
as accurately as possible and calculate the transition scale M c (or rather 
G c {k)). We then calculate the free energy of the cr-model in each phase 
at temperature T and extract the minimal free energy. The deconfinement 
transition occurs at the temperature for which the free energies calculated 
in the ordered and disordered phases of the sigma model coincide. 
The Hclmholtz free energy F of the density matrix p is given by 

F=(H)-T(S), (12) 

where H is the standard Yang-Mills Hamiltonian 



H = / d'> 

with 



l -Ef + l -Bf 



(13) 



5Af(x) ' 

B?(x) = le ijk {djAi(x) - d k A«(x) + 9 r bc A){x)Al{x)} , (14) 



S is the entropy, and T is the temperature. 
Thus 

F = i(Tr(£ 2 g) + Tr(B 2 e)) + T -Trying). (15) 

In the disordered phase, no progress seems possible without restrict- 
ing the arbitrary kernels. Following [4], we adopt the forms cq. (5). For 
small H, we consider only the first non-trivial order in H, that is a term 
of o(HhiH) in the entropy. This term can be written as a product of 
left SU(N) and right SU(N) currents and does, therefore, vanish in the 
disordered, SU(N) L (g> SU(N) R symmetric, phase [4]. The remaining con- 
tribution to the free energy, the average of the Hamiltonian, is evaluated 
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in the mean field approximation [2]. The free energy is minimized for 
M = M c ~ 1.33GeV 

N 2 M 4 

The simplest option to evaluate the free energy in the disordered phase 
is to use perturbation theory. Perturbation theory is certainly appropriate 
for large enough values of M, where the expectation value of the U field is 
close to unity. From numerical studies [10] it is known that the transition 
occurs when the expectation value of U is greater than 0.5. We can thus 
expect perturbation theory to be qualitatively reliable down to the transi- 
tion point. In the leading order perturbation theory approximation to the 
ordered phase of the cr-model, however, minimisation with respect to arbi- 
trary kernels G _1 and H for both high and low modes is possible. Further, 
the analysis can be carried out to all orders in the thermal disorder kernel 
H. 

In this approximation, the U matrices can be parameterized in the stan- 
dard exponential form and expanded in the coupling g 

[7 = cxp|z^ a y| = l + W a y + ... (17) 
Hence at leading order one can take 

U~l, diU-igdi^. (18) 
Thus, the gauge transformations (2) reduce to 

A1 -► A1 - d t <p a (19) 
and the Hamiltonian (13) reduces to 

H=\ [Ef + (e l3k d 3 A a k f] . (20) 

These last two equations describe the theory [/(l)^ 2-1 : in the leading order 
of cr-model perturbation theory, the SU (N) Yang-Mills theory reduces to 
the U(1) N _1 free theory. Moreover, the density matrix eq. (1) becomes 
Gaussian again, because the gauge transformations are linear: 

p[A, A] = J exp | - i [AG^A + (A' - d^p)G- 1 {A l - dip) 

- 2AH(A' - dip)] \ . (21) 
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The theory of TV 2 — 1 U(l) free fields in 3 + 1 dimensions is completely 
tractable; the variational analysis for the U(l) theory (with Gaussian ansatz 
(21)) was discussed in [12]. The free energy in momentum space in terms 
of the arbitrary kernels G _1 and H is 



N 2 



-I 



d 3 p 
(2tt)3 



G _1 (l + GH) + p 2 G(l - GH)~ 



- 4T In 



'GH' 


- In 


V 




[ £ . 




GH 


'V. 



(22) 



where rj = 1 - (1 - {GH) 2 ) 1 / 2 and £ = (1 - (GH) 2 ) 1 / 2 - (1 - GH). The 
kernels which minimize the free energy are 



G- 1 = p 



1 



1-e" 



H = 2p 



e t 



1 - e 



and the minimal value of the free energy at temperature T is 

F = 



T 2 N 2 T 4 



45 



(23) 



(24) 



So we see that the free energy of SU(N) is minimised with M = M c 
in the disordered phase of the cr-model for temperatures from zero up to a 
temperature T c where 



N 2 M* 
' 30tt 2 



which in turn implies 



1/4 



7T 



2 N 2 T? 
1$ 



~ 470MeV . 



(25) 



(26) 



Although the actual value of the transition temperature is considerably 
larger than the lattice estimate it makes more sense to look at dimensionlcss 
quantities. In particular, if we identify 2M C with the mass of the lightest 
glucball [11] (see also the discussion in [1]), we find 



2M r 



1 

2^ 



1/4 



0.18. 



(27) 



This is in excellent agreement with the lattice estimate for SU(3) pure 
gauge theory [8] . 
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4. Conclusions 

In 3+1 dimensional gluodynamics, this variational method gives results 
which on the qualitative level at least, conform with our intuition about 
the structure of the ground state. We find a first order phase transition 
which corresponds to the Polyakov loop acquiring a non-zero average. Al- 
though we have not calculated the string tension directly, the behavior of 
the Polyakov loop is very much indicative that this is indeed the deconfining 
phase transition. The value of the critical temperature (in units of glueball 
mass) we find is in good agreement with lattice results. We also found that 
in the low temperature phase the entropy remains zero all the way up to 
the transition temperature. This is a rather striking result, which has not 
been built into our variational ansatz, but rather emerged as the result of 
the dynamical calculation. 

An important lesson we learned is that the projection of the Gaussian 
trial state onto the gauge invariant Hilbert subspace dictates most, if not all, 
of the important aspects of the non-perturbative physics. It was absolutely 
essential to perform the projection non-perturbatively, fully taking into 
account the contribution of the overlap between gauge rotated Gaussians 
into the variational energy prior to minimization. We have seen that from 
the point of view of the effective c-model the energy is minimized in the 
disordered phase. In other words, the low momentum fluctuations of the 
field U are large, unlike in the perturbative regime, where U is close to a 
unit matrix. From the point of view of the trial wave functional, this means 
that the off-diagonal contributions, coming from the Gaussian WF gauge 
rotated by a slowly varying gauge transformation, are large. It is these 
"off diagonal" contributions to the energy that lowered the energy of the 
best trial state below the perturbative value. In the low temperature phase 
the vanishing of the entropy was also a direct consequence of the effective 
cr-model being in the disordered phase, and thus of the non-perturbative 
nature of the gauge projection. 

This variational method appears to be a good candidate for a useful cal- 
culational scheme for strongly interacting gauge theories. However, many 
outstanding questions remain. Is the best variational state confining? How 
do we calculate the interaction potential between external sources? How 
do we understand better the relation between the variational parameter 
and the glueball masses? Can we extend the Ansatz to include (massless) 
fcrmions? 

Personally, I believe that these results are genuine and that there is 
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enough scope for further development of the approach which warrants con- 
tinuing active investigations. 



Acknowledgments 

This talk is based upon work done in collaboration with Ben Gripaios , 
Alex Kovner and Ian Kogan. 

References 

1. A. Kovner and J. G. Milhano, in M. Shifman, J. Wheater and A. Vainstein 
(eds.): Circumnaviagting Theoretical Physics (World Scientific, Singapore, 
2004), [hep-ph/0406165]. 

2. I. Kogan and A. Kovner, Phys. Rev. D52, 3719 (1995); e-Print Archive: hep- 
th/9408081; 

3. R. Feynman, in Wangerooge 1987, Proceedings, Variational calculations in 
quantum field theory, L. Polley and D. Pottinger, eds (World Scientific, Sin- 
gapore, 1988). 

4. I. I. Kogan, A. Kovner and J. G. Milhano, JHEP 0212, 017 (2002) [arXiv:hep- 
ph/0208053]. 

5. B. M. Gripaios and J. G. Milhano, Phys. Lett. B 564 (2003) 104 [arXiv:hep- 
ph/0302172]. 

6. A. Polyakov, Phys. Lett. B 72 (1978) 477 

7. L. Susskind, Phys. Rev. 20 (1979) 2610 

8. M. J. Teper, arXiv:hep-th/9812187. 

9. W. E. Brown and I. I. Kogan, Int.J.Mod.Phys.A14:799,1999; e-Print Archive: 
hep-th/9705136; 

10. J. B. Kogut, M. Snow and M. Stone Nucl.Phys.B200:211,1982 

11. B. M. Gripaios, Int. J. Mod. Phys. A 18 (2003) 85 [arXiv:hep-ph/0204310]. 

12. B. M. Gripaios, Phys. Rev. D 67, 025023 (2003) [arXiv:hep-th/0211104]. 



